Malaria Mapping Using Transmission Models: Application to Survey Data from Mali

Geographic mapping of the distribution of malaria is complicated by the limitations of the available data. The most widely available data are from prevalence surveys, but these surveys are generally carried out at arbitrary locations and include nonstandardized and overlapping age groups. To achieve comparability between different surveys, the authors propose the use of transmission models, particularly the Garki model, to convert heterogeneous age prevalence data to a common scale of estimated entomological inoculation rates, vectorial capacity, or force of infection. They apply this approach to the analysis of survey data from Mali, collected in 1965-1998, extracted from the Mapping Malaria Risk in Africa database. They use Bayesian geostatistical models to produce smooth maps of estimates of the entomological inoculation rates obtained from the Garki model, allowing for the effect of environmental covariates. Again using the Garki model, they convert kriged entomological inoculation rates values to age-specific malaria prevalence. The approach makes more efficient use of the available data than do previous malaria mapping methods, and it produces highly plausible maps of malaria distributio

Domestic dog demographic structure and dynamics relevant to rabies control planning in urban areas in Africa: the case of Iringa, Tanzania

<p>Background Mass vaccinations of domestic dogs have been shown to effectively control canine rabies and hence human exposure to rabies. Knowledge of dog population demography is essential for planning effective rabies vaccination programmes; however, such information is still rare for African domestic dog populations, particularly so in urban areas. This study describes the demographic structure and population dynamics of a domestic dog population in an urban sub-Saharan African setting. In July to November 2005, we conducted a full household-level census and a cross-sectional dog demography survey in four urban wards of Iringa Municipality, Tanzania. The achievable vaccination coverage was assessed by a two-stage vaccination campaign, and the proportion of feral dogs was estimated by a mark-recapture transect study.</p> <p>Results The estimated size of the domestic dog population in Iringa was six times larger than official town records assumed, however, the proportion of feral dogs was estimated to account for less than 1% of the whole population. An average of 13% of all households owned dogs which equalled a dog:human ratio of 1:14, or 0.31 dogs per household or 334 dogs km-2. Dog female:male ratio was 1:1.4. The average age of the population was 2.2 years, 52% of all individuals were less than one year old. But mortality within the first year was high (72%). Females became fertile at the age of 10 months and reportedly remained fertile up to the age of 11 years. The average number of litters whelped per fertile female per year was 0.6 with an average of 5.5 pups born per litter. The population growth was estimated at 10% y-1.</p> <p>Conclusions Such high birth and death rates result in a rapid replacement of anti-rabies immunised individuals with susceptible ones. This loss in herd immunity needs to be taken into account in the design of rabies control programmes. The very small proportion of truly feral dogs in the population implies that vaccination campaigns aimed at the owned dog population are sufficient to control rabies in urban Iringa, and the same may be valid in other, comparable urban settings.</p&gt

Spatio-temporal modelling of changes in air pollution exposure associated to the COVID-19 lockdown measures across Europe

The lockdown and related measures implemented by many European countries to stop the spread of the SARS-CoV-2 virus (COVID-19) pandemic have altered the economic activities and road transport in many cities. To rigorously evaluate how these measures have affected air quality in Europe, we developed Bayesian spatio-temporal (BST) models that assess changes in the surface nitrogen dioxide (NO2) and fine particulate matter (PM2.5) concentration across the continent. We fitted BST models to measurements of the two pollutants in 2020 using a lockdown indicator covariate, while accounting for the spatial and temporal correlation present in the data. Since other factors, such as weather conditions, local combustion sources and/or land surface characteristics may contribute to the variation of pollutant concentrations, we proposed two model formulations that allowed the differentiation between the variations in pollutant concentrations due to seasonality from the variations associated to the lockdown policies. The first model compares the changes in 2020, with the ones during the same period in the previous five years, by introducing an offset term, which controls for the long-term average concentrations of each pollutant during 2014-2019. The second approach models only the 2020 data, but adjusts for confounding factors. The results indicated that the latter can better capture the lockdown effect. The measures taken to tackle the virus in Europe reduced the average surface concentrations of NO2 and PM2.5 by 29.5% (95% Bayesian credible interval: 28.1%, 30.9%) and 25.9% (23.6%, 28.1%), respectively. To our knowledge, this research is the first to account for the spatio-temporal correlation present in the monitoring data during the pandemic and to assess how it affects estimation of the lockdown effect while accounting for confounding. The proposed methodology improves our understanding of the effect of COVID-19 lockdown policies on the air pollution burden across the continent

Estimating the burden of malaria in Senegal : Bayesian zero-inflated binomial geostatistical modeling of the MIS 2008 data

The Research Center for Human Development in Dakar (CRDH) with the technical assistance of ICF Macro and the National Malaria Control Programme (NMCP) conducted in 2008/2009 the Senegal Malaria Indicator Survey (SMIS), the first nationally representative household survey collecting parasitological data and malaria-related indicators. In this paper, we present spatially explicit parasitaemia risk estimates and number of infected children below 5 years. Geostatistical Zero-Inflated Binomial models (ZIB) were developed to take into account the large number of zero-prevalence survey locations (70%) in the data. Bayesian variable selection methods were incorporated within a geostatistical framework in order to choose the best set of environmental and climatic covariates associated with the parasitaemia risk. Model validation confirmed that the ZIB model had a better predictive ability than the standard Binomial analogue. Markov chain Monte Carlo (MCMC) methods were used for inference. Several insecticide treated nets (ITN) coverage indicators were calculated to assess the effectiveness of interventions. After adjusting for climatic and socio-economic factors, the presence of at least one ITN per every two household members and living in urban areas reduced the odds of parasitaemia by 86% and 81% respectively. Posterior estimates of the ORs related to the wealth index show a decreasing trend with the quintiles. Infection odds appear to be increasing with age. The population-adjusted prevalence ranges from 0.12% in Thille-Boubacar to 13.1% in Dabo. Tambacounda has the highest population-adjusted predicted prevalence (8.08%) whereas the region with the highest estimated number of infected children under the age of 5 years is Kolda (13940). The contemporary map and estimates of malaria burden identify the priority areas for future control interventions and provide baseline information for monitoring and evaluation. Zero-Inflated formulations are more appropriate in modeling sparse geostatistical survey data, expected to arise more frequently as malaria research is focused on eliminatio

Modeling the effect of different drugs and treatment regimen for hookworm on cure and egg reduction rates taking into account diagnostic error

BACKGROUND: Hookworm infections, caused by Ancylostoma duodenale and Necator americanus, are of considerable public health importance. The World Health Organization recommends preventive chemotherapy as the key strategy for morbidity control. Meta-analyses have been conducted to estimate treatment efficacy of available drugs and drug combinations. However, in most studies, the relation between the diagnostic error and infection intensity have not been considered, resulting in an overestimation of cure rates (CRs). METHODOLOGY: A Bayesian model was developed to compare the 'true' CR and egg reduction rate of different treatment regimens for hookworm infections taking into account the error of the recommended Kato-Katz thick smear diagnostic technique. It was fitted to the observed egg count data which was linked to the distribution of worms, considered the day-to-day variation of hookworm egg excretion and estimated the infection intensity-dependent sensitivity. The CR was obtained by defining the prevalence of infection at follow-up as the probability of having at least one fertilized female worm. The model was applied to individual-level egg count data available from 17 treatments and six clinical trials. PRINCIPAL FINDINGS: Taking the diagnostic error into account resulted in considerably lower CRs than previously reported. Overall, of all treatments analyzed, mebendazole administered in six dosages of 100 mg each was the most efficacious treatment with a CR of 88% (95% Bayesian credible interval: 79-95%). Furthermore, diagnostic sensitivity varied with the infection intensity and sampling effort. For an infection intensity of 50 eggs per gram of stool, the sensitivity is close to 60%; for two Kato-Katz thick smears it increased to approximately 76%. CONCLUSIONS/SIGNIFICANCE: Our model-based estimates provide the true efficacy of different treatment regimens against hookworm infection taking into account the diagnostic error of the Kato-Katz method. Estimates of the diagnostic sensitivity for different number of stool samples and thick smears are obtained. To accurately assess efficacy in clinical trials with the Kato-Katz method, at least two stool samples on consecutive days should be collected

Bayesian geostatistical modelling for mapping schistosomiasis transmission

Progress has been made in mapping and predicting the risk of schistosomiasis using Bayesian geostatistical inference. Applications primarily focused on risk profiling of prevalence rather than infection intensity, although the latter is particularly important for morbidity control. In this review, the underlying assumptions used in a study mapping Schistosoma mansoni infection intensity in East Africa are examined. We argue that the assumption of stationarity needs to be relaxed, and that the negative binomial assumption might result in misleading inference because of a high number of excess zeros (individuals without an infection). We developed a Bayesian geostatistical zero-inflated (ZI) regression model that assumes a non-stationary spatial process. Our model is validated with a high-quality georeferenced database from western Côte d'Ivoire, consisting of demographic, environmental, parasitological and socio-economic data. Nearly 40% of the 3818 participating schoolchildren were infected with S. mansoni, and the mean egg count among infected children was 162 eggs per gram of stool (EPG), ranging between 24 and 6768 EPG. Compared to a negative binomial and ZI Poisson and negative binomial models, the Bayesian non-stationary ZI negative binomial model showed a better fit to the data. We conclude that geostatistical ZI models produce more accurate maps of helminth infection intensity than the spatial negative binomial one

Spatial and temporal dynamics of malaria transmission in rural western Kenya

ABSTRACT: BACKGROUND: Understanding the impact of reducing Plasmodium falciparum malaria transmission requires estimates of the relationship between health outcomes and exposure to infectious mosquitoes. However, measures of exposure such as mosquito density and entomological inoculation rate (EIR) are generally aggregated over large areas and time periods, biasing the outcome-exposure relationship. There are few studies examining the extent and drivers of local variation in malaria exposure in endemic areas. METHODS: We describe the spatio-temporal dynamics of malaria transmission intensity measured by mosquito density and EIR in the KEMRI/CDC health and demographic surveillance system using entomological data collected during 2002-2004. Geostatistical zero inflated binomial and negative binomial models were applied to obtain location specific (house) estimates of sporozoite rates and mosquito densities respectively. Model-based predictions were multiplied to estimate the spatial pattern of annual entomological inoculation rate, a measure of the number of infective bites a person receive per unit of time. The models included environmental and climatic predictors extracted from satellite data, harmonic seasonal trends and parameters describing space-time correlation. RESULTS: Anopheles gambiae s.l was the main vector species accounting for 86% (n=2309) of the total collected mosquitoes with the remainder being Anopheles funestus. Sixty eight percent (757/1110) of the surveyed houses had no mosquitoes. Distance to water bodies, vegetation and day temperature were significantly associated with mosquito density. Overall annual point estimates of EIR were 6.7, 9.3 and 9.6 infectious bites per annum for 2002, 2003 and 2004 respectively. Monthly mosquito density and EIR varied over the study period peaking in May during the wet season. The predicted and observed densities and EIR showed a strong seasonal and spatial pattern over the study area. CONCLUSIONS: Spatio-temporal maps of malaria transmission intensity obtained in this study are not only useful in understanding variability in malaria epidemiology over small areas but also provides a high resolution exposure surface that can be used to analyse the impact of malaria exposure on mortalit

Constructing a malaria-related health service readiness index and assessing its association with child malaria mortality: an analysis of the Burkina Faso 2014 SARA data

Background: The Service Availability and Readiness Assessment surveys generate data on the readiness of health facility services. We constructed a readiness index related to malaria services and determined the association between health facility malaria readiness and malaria mortality in children under the age of 5 years in Burkina Faso. Methods: Data on inpatients visits and malaria-related deaths in under 5-year-old children were extracted from the national Health Management Information System in Burkina Faso. Bayesian geostatistical models with variable selection were fitted to malaria mortality data. The most important facility readiness indicators related to general and malaria-specific services were determined. Multiple correspondence analysis (MCA) was employed to construct a composite facility readiness score based on multiple factorial axes. The analysis was carried out separately for 112 medical centres and 546 peripheral health centres. Results: Malaria mortality rate in medical centres was 4.8 times higher than that of peripheral health centres (3.5% vs. 0.7%, p < 0.0001). Essential medicines was the domain with the lowest readiness (only 0.1% of medical centres and 0% of peripheral health centres had the whole set of tracer items of essential medicines). Basic equipment readiness was the highest. The composite readiness score explained 30 and 53% of the original set of items for medical centres and peripheral health centres, respectively. Mortality rate ratio (MRR) was by 59% (MRR = 0.41, 95% Bayesian credible interval: 0.19-0.91) lower in the high readiness group of peripheral health centres, compared to the low readiness group. Medical centres readiness was not related to malaria mortality. The geographical distribution of malaria mortality rate indicate that regions with health facilities with high readiness show lower mortality rates

Spatial Patterns of Infant Mortality in Mali: The Effect of Malaria Endemicity

A spatial analysis was carried out to identify factors related to geographic differences in infant mortality risk in Mali by linking data from two spatially structured databases: the Demographic and Health Surveys of 1995-1996 and the Mapping Malaria Risk in Africa database for Mali. Socioeconomic factors measured directly at the individual level and site-specific malaria prevalence predicted for the Demographic and Health Surveys' locations by a spatial model fitted to the Mapping Malaria Risk in Africa database were examined as possible risk factors. The analysis was carried out by fitting a Bayesian hierarchical geostatistical logistic model to infant mortality risk, by Markov chain Monte Carlo simulation. It confirmed that mother's education, birth order and interval, infant's sex, residence, and mother's age at infant's birth had a strong impact on infant mortality risk in Mali. The residual spatial pattern of infant mortality showed a clear relation to well-known foci of malaria transmission, especially the inland delta of the Niger River. No effect of estimated parasite prevalence could be demonstrated. Possible explanations include confounding by unmeasured covariates and sparsity of the source malaria data. Spatial statistical models of malaria prevalence are useful for indicating approximate levels of endemicity over wide areas and, hence, for guiding intervention strategies. However, at points very remote from those sampled, it is important to consider prediction erro

The economic impact of schistosomiasis

BACKGROUND: The economic impact of schistosomiasis and the underlying tradeoffs between water resources development and public health concerns have yet to be quantified. Schistosomiasis exerts large health, social and financial burdens on infected individuals and households. While irrigation schemes are one of the most important policy responses designed to reduce poverty, particularly in sub-Saharan Africa, they facilitate the propagation of schistosomiasis and other diseases. METHODS: We estimate the economic impact of schistosomiasis in Burkina Faso via its effect on agricultural production. We create an original dataset that combines detailed household and agricultural surveys with high-resolution geo-statistical disease maps. We develop new methods that use the densities of the intermediate host snails of schistosomiasis as instrumental variables together with panel, spatial and machine learning techniques. RESULTS: We estimate that the elimination of schistosomiasis in Burkina Faso would increase average crop yields by around 7%, rising to 32% for high infection clusters. Keeping schistosomiasis unchecked, in turn, would correspond to a loss of gross domestic product of approximately 0.8%. We identify the disease burden as a shock to the agricultural productivity of farmers. The poorest households engaged in subsistence agriculture bear a far heavier disease burden than their wealthier counterparts, experiencing an average yield loss due to schistosomiasis of between 32 and 45%. We show that the returns to water resources development are substantially reduced once its health effects are taken into account: villages in proximity of large-scale dams suffer an average yield loss of around 20%, and this burden decreases as distance between dams and villages increases. CONCLUSIONS: This study provides a rigorous estimation of how schistosomiasis affects agricultural production and how it is both a driver and a consequence of poverty. It further quantifies the tradeoff between the economics of water infrastructures and their impact on public health. Although we focus on Burkina Faso, our approach can be applied to any country in which schistosomiasis is endemic